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1. 


Introduction 


The thermal transport properties of three dimensional carbon-fibre 
composites (C.F.C.C.) have been measured over the temperature range 
300-3000K. 

The thermal diffusivity (a) has been measured over this temperature 

(1 2 ) 

range using the laser pulse technique. * Data thus obtained has been 
converted to thermal conductivity (X) using the standard expression 
X = a p Cp 

where p is density and Cp is specific heat. 

The structure of these materials is such that it is necessary to 

make measurements along two axes in order to characterise the properties. 

Although the experimental data is, in its own right extremely 
valuable, an additional aim of this investigation has been an attempt 
to model the C.F.C.C. thermal conductivity in terms of the properties of 
the constituents. Three different 3-D composites were supplied by the 
A.F.M.L., Wright-Patterson A.F.B. These were 

a) Material A 3-D orthogonal C.F.C.C. using Fibre 

"F" material. 

b) Material B 3-D C.F.C.C. using Fibre 

"F" material using on 8-harness satin weave. 

/ 

c) Material C- 3-D orthogonal C.F.C. using Fibre "A” High 

Modulus(H.M.) fibre material. 

Samples of matrix material and 1-dimensional composites to assist modelling 
were also obtained and measured over the same temperature range. 

Although the finished 3-D C.F.C.C. consists only of a porous 
body containing carbon fibre yarns and graphitised matrix material there 
are many variables within processing so an accurate description of the 
finished composite may be much more complex. A fairly simple modelling 
approach has been adapted since accurately quantifiable information on 
constituent properties, from 1-D composites and matrix material is not 


available. 







Data for composites, A and B (fibre F) has been detailed in 
the progress report for the period ending 30/9/79. Consequently this 
report will concern itself primarily with data for material "C" and 
comparable 1-D composites and an evaluation of comparisons with data 
for materials "A" and "B". 







2. Characterisation of Materials 

Materials "A" and "C" are 3-D orthogonal C.F.C.C. in which the carbon 
fibres are formed into yarns. In the case of composite "A” this is 
known as a 2-2-3 composite comprising two yarns of 1440 filament yarns in the 
X and Y direction and three yarns of 1440 filaments forming the reinforcement 
in the Z direction. A unit cell of this material is shown in figure 1. 
Composite "C" is similar but in this case is a 2-2-1 composite. The X and 
Y reinforcement is composed of two yarns of 1000 filaments/site whereas 
the Z axis reinforcement is composed of one 3000 filament yarn. Material 
B is also a 3-D composite, but in this case the yarns are not orthogonal. 

The X and Y yarns are woven into an 8-harness satin weave to form a rein¬ 
forcement plane which is pierced by the Z-yarns. This type of 
construction is shown in figure 2. 

The manufacture of a typical composite begins with the construction 
of a fibre preform. The second stage is to coat the preform with 
pyrolytic carbon. This infiltration of the preform is done using low 
pressure, isothermal deposition of vapour phase carbon (C.V.D) prior to 
densification of the preform by pitch impregnation, and high pressure 
pyrolysis. The C.V.D. stage is instituted to stiffen the preform and 
prevent deformation during high pressure pyrolysis. This impregnation 
process may be typically repeated 5 times. The composite is finally 
graphitised at 3023K (2750°C). The objective of this processing 
is to impregnate the fibre bundles and fill crossover pockets and voids 
with graphitic matrix material. 

The following additional specimens were provided by the A.F.M.L. 

Wright Patterson A.F.B. 

1) Samples of Fibre "A” and Fibre "F" material used 
in the composites. 

2) A sample comprised only of graphitised matrix material. 

3) A 1-D composite made using fibre "F" material that 
had not been subject to the C.V.D. process prior to 
densification and graphitisation. 










4) 


A 1-D composite made using Fibre "A" material that had been 


C.V.D.'d and graphitised to 2573K (2300°C). 

5) A 1-D composite made using Fibre "A" that had been 
graphitised to 2973K (2700°C) but not C.V.D.’d. 








Table 1. Materials Specification 





























3. Measurement Technique 

The U.M.I.S.T. apparatus (figure 3) has recently been described 

( 3 ) 

in detail so will only be briefly described here. The front face 
of a disc shaped sample is heated instantaneously 10 3 s) using a 
100J ruby laser. The specimen is heated to the measurement temperature 
inside a graphite susceptor heated using an induction coil. The 
subsequent temperature transient on the rear face of the sample is 
recorded using an infra-red sensor. Thermal diffusivity is 
calculated from 


a 




( 1 ) 


where L is the sample length. 

t^ is the time taken to attain half the maximum temperature rise 
w 

—^ 5 * is a constant, which in the absence of heat losses = 0.139. 
The amplified signal is digitally logged by a microcomputer which 
samples data at 0.001s intervals. 

Heat losses are inevitable at high temperatures and necessitate 

w (4) 

that —£ decreases. This has been analysed by Cowan . To minimise 

this correction and avoid overly large corrections for finite pulse time 

affects it is desirable to select a sample length such that 

0.025s < t, < 0.2s.Suitable programming enables t •, to'be computed to 
5 2 

precision better than 10 4 s (including correction for finite pulse 

time where necessary). The ratio of.transient amplitudes AT(10t,) 

i 

/AT(ti) or AT(5t i) /AT(t-i) is then determined and —r calculated using 
the analysis outlined by Cowan. Finally the sample length is corrected 
for thermal expansion and the diffusivity computed from these three 


parameters. 






4. 


Specimens G.35mm diameter and of length chosen to fit the experimental 
requirements listed in the preceding section were cored from each 
material billet provided. Samples were cored from the Z and X directions 
of the 3-D C.F.C.C. materials since they are symmetrical in the X and 
Y directions. Specimens perpendicular and transverse were machined from 
the 1-D composites. For the matrix graphite the isotropy of the material 
was verified using samples cored from two orthogonal directions. Sample 
data is presented in Table 2. 

Density measurements were obtained using a liquid densitometer 
and reveal only minor variations between samples. Bulk densities are 
in general 2-3% lower which indicates a small amount of open porosity. 
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Composite 

Type 


3D CFCC "A" 


3D CFCC "B" 
Fibre "F" 


1-D CFCC 
Fibre "F" 
(2750°C) 


Matrix 

3-D CFCC "C" 


1-D C FCC 
Fibre "A" 

CVD (2300°C) 

1-D C FCC 
Fibre "A " 

No CVD 
(2750°C) 


Table 2. Specimen Details 


Orientation 

Sample 

No 

Length 

cm 

Density 
gmcm 3 






X axis 

* 

1 

0.242 

1.926 



. 2 

0.242 

1.920 

< Z axis 


r 1 

0.399 

1.916 



2 

0.398 

1.887 






I 

3 

0.298 

1.919 

' i 

l i 

4 

0.301 

1.925 

1 

X axis , 


0.439 


__1 

L 2 

0.439 



Z .ixls 


Parallel to 
fibre axis 
Transverse to 
fibre axis 


1 

2 

3 

1 

3 

1 


X axis 


1 

2 

1 

2 

3 


0.338 

0.336 

0.338 

0.454 

0.458 

0.1995 

0.152 

0.1992 

0.297 

0.3005 

0.3068 

0.257 

0.3105 


Z axis 

Parallel to fibre 
axis 


Transverse to 
fibre axis 

Parallel to 
fibre axis 


Transverse to 
fibre axis 


i 1 

0.3668 

1.880 


0.3601 

1.888 

X 

0.4039 

1.622 

2 

0.353 

1.626 

1 

0.2065 

1.629 

2 

0.2135 

1.624 

1 

0.4086 

1.665 

2 

0.3561 

1.668 

1 

0.2128 

1.701 

2 

0.2012 

1.665 


1.936 

1.931 

1.928 

1.910 

1.905 

1.900 

1.897 

1.910 

1.369 

1.389 


1.890 
1.876 
1.886 




















5. Experimental Procedur e 

Diffusivity measurements were carried out on all the materials 
listed in section 2 (except for the Fibre "A" 1-D composite graphitised 
to 2573K) over the temperature range 300-3000K. It is necessary to make 
these measurements in four distinct stages as follows. 

a) 350-750K. Transients were recorded using a Hg^ Cd^_ x Te 
detector (cut off wavelength 5.5 gm) and a synchronous (detector) 
lock-in amplifier. This was an experimental system replacing the 
earlier InSb detector (used for fibre "F" composites ) and was only 
useable down to approximately 350K at acceptable noise levels. 

b) 500-1500K. Transients were detected using PbS detector (cut 
off wavelength 3.0 pm) and the measurement temperature recorded using 
a Type K thermocouple. 

c) 1200-2300K. As b) above but temperatures recorded using an 
optical pyrometer. 

d) 2000-3000K. A vacuum (< 10 4 Torr) was sufficient to protect the 

samples over the preceding temperature ranges. Above 2300K the vapour 
pressure of graphite increases so it is necessary to suppress evaporation 
This was done by completing measurements in an atmosphere of helium 

at 20 p.s.i. pressure. 

The above temperature ranges all overlap. This is to confirm that 
these changes in measurement technique did not yield different 
diffusivity values. 

For modelling purposes it is necessary to convert the measured 
thermal diffusivity to thermal conductivity. This requires accurate 
data for specific heat and thermal expansion. The specific heat data 

t C \ 

used was supplied by A.F.M.L. synthesised from measurements from 

( 8 ) 

... (7) and a C.F.C.C. composite "A" from 

350-1000K on P0C0 graphite r 

1500-3000K. This data is fitted to a 5 th order polynomial using a 
least squares function (figure 3). 




Cp = -5.444xlO _1 + 5.5076x10" 3 T - 4.9454xl0 _6 T 2 

+2.3389x 10 _9 T 3 + 5.5749xl0 _1 3 T l( - 5.3241xl0 _17 T 5 

( 0 ) 

Thermal expansion data supplied by A.F.M.L. has been used to 
correct for density changes with temperature. 

Specimens for microstructural investigation were prepared by 
grinding and polishing down to 1/4 pm. These were then etched for 
1-5 minutes in a hot solution of potassium dichromate in phosphoric 
acid. In order to improve resolution in the electron microscope, 
the etched specimens were then lightly sputtered with gold. 
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Results 


At least two specimens were measured for each type of specimen. 
Conductivity values derived from the diffusivity data were fitted 
with least squares functions. A number of different functions 
were tried and the best fit over the whole temperature range was 
usually a fourth order polynomial. Although data for each specimen type 
is presented briefly in section 6(A) - 6(G) the coeficients and root mean 
square errors for the best fit functions are listed in Table 3. 

(A) 3-D Composite 'A' 

(i) X Axis 

The diffusivity of two X axis samples was measured over the complete 

temperature range. The results are shown in Figure 4. They provide an 

excellent agreement with AFML flash diffusivity data, also shown, which 

was obtained from a CFCC similar to ’A'. 

Figure 5 shows the X axis conductivity together with directly 

measured data obtained from different billets of 'A' type material by 

( 6 ) 

the comparative rod method' . The direct method conductivity is 
higher at low temperatures, and this may simply be a product of inter¬ 
billet variations. However Lee and Taylor^^ and Minges^ 10 ^ have 
observed similar discrepancies between direct and diffusivity-derived 
conductivity. 

There is no systematic difference in the properties of the two 
samples and this is shown more clearly in Figure 6 which is a plot of 
% deviation from the least squares function. 

(ii) Z Axis 

Figure 7 shows the diffusivity of four samples, again compared 
with diffusivity data from a material similar to composite ’A'^ 6 ^. 

As in the previous case, the agreement of the diffusivity data is good. 

Figure 8 shows the resultant conductivity values with comparative 
rod data for comparison^. As before, the latter data is higher. Billet 

(9) 

variations or, as Lee and Taylor imply, differences in measurement 














































































method may be responsible. In this case, however, it seems appropriate 
to question the validity of the low temperature specific heat data. Over 
a temperature range of 350-1000K this relies entirely on the specific 
heat of POCO graphite, whereas there is evidence that the specific heat 
of fibre 'F' alone is 30% higher at 300K . The law of mixtures would 

suggest that the composite value should lie between these two limits. 

Figure 9 shows the data deviation plot for the least squares function 
given in Table 3. A systematic relationship between deviation and sample 
can be seen which may be partly, but not completely, explained by 
differences in sample density (see Table 2). 

(B) 3-D Composite 'B' 

(i) X Axis 

Diffusivity data and derived conductivity values of two X axis samples 
are shown in Figures 10 and 11 respectively. No comparative diffusivity 
data is available. However, comparative rod conductivity data v , 
from a different billet, shows excellent agreement. 

Figure 12 shows the deviation from the least squares fit. 

(ii) Z Axis 

Diffusivity and conductivity results from three samples are shown 
in Figures 13 and 14 respectively. Comparative rod data^ & \ again 
from a different billet, is higher by a factor of 20-25% over most of the 
temperature range. 

Figure 15 shows that deviation from the least squares fit is 
broadly sample dependent. Differences in sample densities are insignificant. 

(C) 3-D Composite ’C 1 
(i) X-Axis 

Diffusivity data and derived conductivity values for three X-axis 


samples are shown in figures 16 and 17. No comparative data is available. 
There was some indication of difference, between sample 1 and samples 
2 and 3. This data is very similar to that obtained for composite 'A'. 








Figure 18 shows the deviation from the least squares fit given 
in Table 3. 

(ii) Z Axis 

Diffusivity data and derived conductivity values for two 
Z-axis samples are shown in figures 19 and 20 respectively. There is 
no detectable sample difference and again similarity to the data obtained 
for composite 'A'. 

Figure 21 shows the deviation from the least squares fit given 
in Table 3. 

(D) Matrix Material 

Diffusivity and conductivity results from two bulk matrix samples 

are shown in Figures 22 and 23 respectively. 

The samples were cut orthogonally to determine the degree of 

anostropy, if any. The data indicate an essentially isotropic material. 

The high porosity of the bulk matrix (^40%) contrasts with the high 

density of the composites and, by inference, with the low porosity of the 

matrix material within the composite. Evidence suggests in general 

that the density of carbon composite phases in situ is higher than when 

processed individually^'*"'*'^. However, the matrix 'crossover’ pockets (see 

Figure 1) are known to have porosity values of between 10-50%^ 6 \ 

The matrix conductivity is very similar to that of P0C0 graphite 

type AXM-5Q1, as well documented reference material^ . With an average 

density of 1.75 grms/cm 3 , the P0C0 poro'sity is considerably lower, 

however. Whilst radiative and gaseous conduction modes may be.come important 

( 12 ) 

in high porosity, low conductivity materials no enhancement of the matrix 

thermal properties was observed in a helium atmosphere both at low and high 
temperatures. Similarly, such expressions as are available to quantify 
raidative conduction in porous materials suggest that this too would be 


insignificant in the present case. 


(E) Fibre "F" 1-D Composite. No CVD -2700 C Graphitisation 


(i) 


Axis Parallel with Fibres 

Diffusivity and conductivity results from two parallel (//) axis 






samples are shown in Figures 24 and 25. No comparative data is available. 

The 1-D composite was unique in that its properties change 
significantly after heating to temperatures > 2300K. Spot measurements, 
at temperatures of 600-700K, were made after the principal measurements 
had been completed. These revealed that the diffusivity of both samples 
was reduced to 50% or less of its original value. In addition there 
was a permanent change in simple length and diameter of + 1-1.5% respectively. 

Since there are no orthogonal yarn influences, it is expected 
that the 1-D composite thermal expansion be higher than that of the 
3-D materials' . However, no permanent offsets have been reported 
before. Significantly, it is now known that the processing of the 

1- D material differs from that of the other composites'' reported here. 

Typically, the high temperature (2000-3000K) measurement runs lasted 

2- 3 hours. Oxygen contamination has been ruled out since other sample types 
were not effected. 

Subsequent liquid densitometer measurements showed a density increase 
of about 1%. Compared with the nominal decrease in bulk density, this 
indicates a decrease in closed porosity and micrstructural changes. 

The high conductivity of the 1-D parallel axis is another surprising 
feature. The room temperature value of approximately 3.0 W/cmK must be 
compared with available data on Fibre ’F’ conductivity which indicates a 
value of = 0.6 W/cmK along the fibre axis^*^’^^. This large difference 
in the properties of the 1-D composite .and fibres has important 
implications and will be discussed more fully when composite modelling 
techniques are examined. 

Doubts have been raised concerning the efficacy of diffusivity 
measurements on highly orientated composites where the integrity of 
conduction paths is preserved. A good example of this is the parallel 
axis of the 1-D composite. Arguably, an effective composite diffusivity 
does not exist for such materials in which event the rear face temperature 
analysis (see Equation 2) will give meaningless results., except in 


two extreme cases. 
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These are: 

u 2 
a 2 

where and are the diffusivities of the phases. 

(9) 

Case 2) has been used by Lee and Taylor to measure the diffusivity 
of carbon fibres. 

The rear face temperature analysis enables the calculation of 
diffusivity with any value of x (0<x<l) in t^provided the appropriate value 
of w/tt 2 is used in Equation 2. Calculation of diffusivity for different 
values of x gives a check on how closely the observed transient conforms 
to its theoretical form and thus may provide information as to the nature 
of the sample. 

The rear face temperature transients of several 1-D composite samples 
were recorded over a temperature range where heat losses would not distort 
the transient. The difference in diffusivity values calculated from t^ 
values in the range 0.2 < t^ < 0.8 was less than ±4%. It may be concluded 
therefore that the parallel axis does have a meaningful effective 
diffusivity. 

(ii) Axis Transverse to Fibres 

Diffusivity and conductivity results for the transverse fibre (TF) 
axis from three samples are shown in Figures 26 and 27. 

The samples transverse to the fibre axis also showed irreversible 
dimensional changes after heating to 2400K. Increases in sample length 
were higher being in the range 2.6-4.9%. Fibre thermal expansion is 
known to be higher in the transverse direction^’^. 

(F) 1-D Composite Fibre "A". No CVD - Graphitised at 2750°C 

(i) Axis parallel to fibres 

Measured thermal diffusivity and derived thermal conductivity 
data for two samples are shown in figures 28 and 29. The diffusivity 
is some 30% lower than the data obtained for the fibre "F" - 1-D composite. 

The density is also significantly lower (1.7 gm cm 3 compared to 1.89 gm cm 3 ) 


1) a. 


>> 
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and the derived conductivity values lie some 37% lower. 

(ii) Axis transverse to fibres 

In similar fashion the thermal diffusivity data obtained for this 
direction (figure 30) are also some 30% lower and the thermal conductivity 
values (figure 31) are likewise 37% lower. 

(G) 1-D Composite Fibre "A 11 + CVD Graphitised at 2300°C 

(i) Axis parallel to fibres 

Thermal diffusivity data for 2 samples are shown in figure 32. 

These data are much lower than those for the 1-D composite of Fibre "A" 
without CVD by a scaling factor of 0.3-0.5. The influence of CVD is almost 
impossible to ascertain since the difference is almost certainly due 
primarily to the difference in graphitisation temperature. An attempt was 
made to graphitic one sample by heating to 3000K. However the diffusivity 
decreased and subsequent examination revealed cracks between fibres 
and between fibres and matrix. 

(ii) Axis transverse to fibres 

Thermal diffusivity data are presented in figure 33 obtained for 
this direction showed values exhibiting a different temperature dependence 
to the values obtained for fibre "A" with no CVD. At T > 1000K values 
are similar but at 300K values are only 70% of the non CVD'd samples. 
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7. Micron lructur.il Invest if,at ion . 

Figures 34 and 35 show scanning electron microscope (C.E.M.) views 
across the fibre axis of two 1-D composite samples, made from fibre "F" 
material. The latter is sample Z after measurements up to 2900K whereas 
the former lias undergone no heating. There were no obvious microstructural 
changes which could explain the large transverse offsett. Fibre/matrix 
interfaces appeared very similar with no evidence of increased separation, 
as higher magnification shows (figures 36 and 37). 

Figure 38 shows a general view across the X axis of an unheated 
specimen of C.F.C.C. "A". The fibre yarns running from bottom left to 
top right are Y axis reinforcement, the other fibre yarns are Z axis. Pores 
of up to lOOwm were visible in the matrix "cross over" pockets. 

Broken yarn/yarn interfaces were a common feature of the 3-D composites 
as illustrated in figure 39 whereas yarn/matrix interfaces were generally 
intact. The X axis reinforcement is composed of two yarns but there was 
little evidence of yarn splitting in the unheated composite. 

Figure 40 shows greater separation between fibres and matrix than 
was visible in the 1-D composite. In some regions the volume between 
adjacent fibres was completely filled with material of an apparently different 

form to the rest of the matrix (see figure 41). This quite possibly is CVD 

+ . .( 6 ) 

material 

A general view across the Z axis of a C.F.C.C. "A" specimen that 
had been measured up to 2900K is shown in figure 42. Splitting within the 
Z axis yarns could be seen and the pattern of broken yarn/yarn interfaces 
repeated but more clearly marked than in the unheated specimen (figure 43). 
However spot diffusivity measurements taken during cooling had revealed no 
measurable change in thermal properties. No evidence could be seen of 
changes in the fibre/matrix interface (figure 44). 

A general view across the X axis of C.F.C.C . "B" is shown in figure 


45, illustrating the different method of fabricating this composite. The 
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X-reinforcements were not easy to distinguish individually because of 
the apparent absence of matrix "cross over" pockets. The area between 
the Z yarns appeared very homogeneous a feature further illustrate d 
in figure 46. Presumably this is a result of using a weave instead of 
orthogonal fibres held in place before CVD by a jig^\ The X and Y 
reinforcements are therefore able to expand into the cross over pockets, 
otherwise filled by the matrix. The vertical cracks, both major and minor 
were regularly spaced and probably mark the boundary between adjacent 
X-Y reinforcement planes which are nominally 0.254mm (254pm thick). 

Broken Z yarn/X-Y interfaces were quite prominent and the fibre/matrix 
microstructure was very similar to that of C.F.C.C. "A". Figure 47 shows 
a general view across the Z axis illustrating the fine weave pierced fabric 
construction. Broken yarn/yarn interfaces are clearly visible. 

However when we examine the composites nominally made from fibre "A" 
then discrepancies become evident. Some very good micrographs were obtained 
of the two 1-D fibre "A" composites which clearly show a circular fibre 
shape. In figure 48 for the 1-D composite that has been CVD’d can be distinctly 
seen a number of spherical fibres of about 6pm diameter. At higher 
magnification in a more heavily etched region of the specimen (figure 49) can 
be seen two fibres which have clearly visible around them a sheath, 
presumably of CVD graphite. In figures 50, 51 and 52 at magnifications of 
x 2,000, x 5,000 and 10,000 are three more regions of the CVD'd 1-D 
composite. Again the graphite sheath around each individual fibre is 
evident but more significantly the structure of the in-fill graphite, presumably 
matrix appears to orient itself transverse (TOG) with respect to the sheath 
of CVD graphite around each fibre which appears to by parallel oriented graphite 
(P.O.G.) 

By way of contrast the non CVD unidirectional composite has a markedly 
different microstructure. In figure 53 is shown a micrograph at x 2,000 
which may be compared with figure 50. Again whilst the cylindrical nature of 
the fibres is clearly evident there is a relative lack of detail regarding 




the microstructure of the matrix. Progression to higher magnification 


figure 54 (x 5,000) and figure 55 (x 10,000) which are to be compared with ' 

j 

figures 51 and 52, reveal little more detail. At x 20,000 magnification 
(figure 5G) some detail is revealed with slight evidence that matrix j 

graphite does form a sheath around the fibres. However no further detail ' 

can be drawn regarding the nature of the microstructure. When a 
longitudinal section of the non CVD'd 1-D fibre "A" composite is examined 
a little more detail is apparent. In figure 57 there is clear evidence that 
the matrix graphite surrounding the fibre is transversely oriented. 

If we now consider the 3-D composite supplied as being made of 
fibre "A" material then a disconcerting difference becomes apparent. In 
figures 58 and 59 are electron micrographs of the Z axis direction and X 
axis direction respectively. Higher magnifications however clearly show 
that the fibres are crenulated (figure 60). This specimen has been heavily 
etched to show the relief. Even on more lightly etched samples figures 
61 (x 5,000) and 62 ( x 10,000) the fibre shape can clearly be seen. In 

these two specimens can also be seen some outcrops of matrix graphite 
which has obvious structural features of similarity with the CVD'd 1-D 
composite (figure 50). 

Figures 63 ( x 200) and 64 ( x 1,000) exhibit evidence of cracks 
at yarn/yarn interfaces whereas there is a better continuity between the 
interfaces between matrix and yarns. There is evidence of substantial 
porosity in the matrix graphite in figure 64. 

It is clear that the 3-D composite is not made from fibre "A" 
unless some extremely unusual etching effect is evident in the 3-D composite 
that is not apparent in the 1-D composite. We must therefore conclude 
that the wrong composite has been supplied and moreover is most probably made 
from fibre "F" material. This does not appear to be composite "A", although 
the diffusivity results are very similar. Analysis of figures 58 and 59 suggests 
that the unit cell dimensions in the Z and X directions are 0.685mm and 

590mm respectively which is substantially less than the unit cell dimensions 

A 




of composite "A." (0.8't x 0.76 mm). This is a scaling factor of 0.8 
and it is extremely unlikely that the electron microscope calibration 
could be in error by 20%. 




8. Modelling of C.F.C.C. Thermal Conducitivty 

The intended approach of determing the thermal transport 

properties of the individual fibre and matrix phases and combining this 

data in the appropriate 3-D geometry is clearly invalidated by the lack 

of suitable specimens from which 1-D composites and matrix graphite 

experimental data may be analysed, and by the fact that composite "C" does 
not appear to be made from fibre "A" material. 

The thermal conductivity of a composite will be a function of 

the thermal conductivities of the matrix material and of the fibres. 

These will be affected by processing variables; matrix graphite by 
the presence of porosity and the fibres by preform stiffening by carbon 
vapour deposition. To that extent none of the back-up materials; the 
three 1-D composites and the matrix graphite, can be said to be represent¬ 
ative of the 3-D composites. 

The bulk density of the matrix graphite is only 1.36 gmem 3 

whereas the bu 1 k densities of the composities are much higher. This 

inplies that some scaling factor is required to predict the conductivity of 

the matrix material of the composites. Again whilst composites A and B and 

the 1-D fibre "F" composite have comparable densities, that for composite 

C and the two fibre "A" 1-D composites differ significantly. Additionally 
the fibre F and fibre A composites graphitised to 2750°C did not have 

CVD treatments and cannot be consider typical of the 3-D composites 

whereas the only 1-D composite that did have CVD treatment was only 

graphitised to 2300°C. Hence it is not possible to take into account 
the influence of processing upon the properties of the constituent phases. 

This may be qualitatively summarised as follows: 

1) Constituent densities tend to be significantly higher within the 

CFCC Increases of 20% and more have been reported in both matrix 

and fibres^^. In the case of the fibres, some of this density increase 

(14) 

is due to the filling of pores in the fibre bundles 

2) If used, the CVD process produces a sheath-like coating around the 
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individual fibres (see Figure 31). The form of the coating varies with 


fibre type and in the case of fibre 'F' is isotropic or 


transverse . Where fibres are closely packed (separation less than l]> 


fibre diameters) the inter-fibre space may be entirely filled with CVD 

_ . ,U7) 

material 


3) The elastic modulus, thermal conductivity and bulk density of CVD 
CFCC are lower than when CVD is absent. The CVD process tends to seal 


off the fine porosity of the fibres, preventing infiltration by the matrix 
(14,15,16) 


4) The matrix crystallites in the vicinity of fibres with isotropic CVD 

(14) 

coatings are transversely oriented (TOG) w.r.t. the fibre axes 

In other regions e.g. cross-over pockets and v/here spacing between fibres 

is large, the matrix is isotropic. 

5) Matrix crystallites in the vicinity of non CVD'd 'F' fibres are 

(’ 7) 

parallel oriented (POG) with the fibre axes 

It is clear that the properties of CFCC materials result from 
a complex interaction of the individual phases. In the case of non 
CVD'd yarns, the matrix must not only increase fibre conductivity by 
infill of fibre defects, but also with its POG structure contribute 
very significantly to axial conductivity. 

In the case of CVD'd material, it is more likely that CVD infill 
is the major influence on axial conductivity with that of the matrix less 
prominent. In all examples however the matrix properties themselves will 
be influenced by FVF since this is a factor in determining the 
proportions of POG/TOG and isotropic matrix. A schematic model of composite 
A is reproduced in figure 65. 
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8.1 Matrix Graphite 

The large difference in density between the as supplied matrix 
graphite and the fibre composites raises doubts as to whether the 
measured values are appropriate for use in composite modelling. 

It has been conclusively shown that the conductivity of graphite is 
dominated by conduction along layer planes and that the thermal 
conductivity of polycrystalline graphite may be represented by ' ’ 


i = £6 r = » r (2 > 

a a 

where X^ is the thermal conductivity along layer planes and the 

parameters e and 8 represent a tortuosity factor and a porosity factor 

( 18 ) 

respectively. Taylor et al have analysed a number of graphites using 
this form of analysis and consider that three scattering processes 
contribute to the thermal resistivity. 



where X^, X^, and X^ represent the contribution to thermal conductivity 
due to Umklepp scattering, grain boundary and isotope scattering. 

Taylor et al have shown that that isotope scattering will contribute 
some 2-3% to the total thermal resistance so, in view of the assumptions 
to be made, this will be neglected. An empirical relationship for the mean 

(19) 

free path has been derived by Taylor 

X u = 8.75xl0" 7 exp 1130 /T * (4) 


from which the thermal conductivity X u = 1/3 CVL> u may be calculated. 
(Values are listed up to 1000K by Taylor). Explicit equations for 
boundary scattering have been derived by Kelly 


K l + 


K 2 + 


(5) 




T 


+ 


T 
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where the symbols No,V ,V ,h, 6 ,c ,0.,0 and W 3 are as used previously 

L» X cl la T 

by Kelly. 

At high temperatures it is reasonable to assume that UmkJapp 
scattering will dominate and that boundary scattering will contribute 
< 10% to the total thermal conductivity. Data for matrix graphite 
can be fitted by a curve for which a = 9.0 and = 2,750 °A. 
to an accuracy of + 8 % over all the curve. This value of crystallite 
size is slightly higher than those derived by Taylor for a range of 
polycrstalline graphites but similar to a value determined by Kelly 
and Gilchrist^ 21 ^ for pyrolytic graphite graphitised at 2750°C. The 
value of 9.0 for the porosity/tortuosity faccor is higher than those 
•noted by Taylor; however the low density (40% porosity) of the bulk 
matrix graphite is probably a prime contributory factor. 

This matrix graphite with a bulk density of 1.36 gm cm 3 has to be 
compared with bulk densities in the range 1.88-1.92 gm cm -3 for 
the two 1-dimensional fibre "A" composites. Hence to use the raw 

derived thermal conductivity data is probably unrealistic. A fully 

dense isotropic graphite would have zero porosity factor and a 

tortuosity factor of /2~ . A line of interpolation between a = 1.41 

at zero porosity and a = 9.0 at 40% porosity should however permit an 

evaluation of a for any value of density. Although there is 'no 

justification for this,it is perhaps a reasonable assumption in view of the 

complex interaction of variables and lack of any available data. 
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8 . 2 . 1-F) Composites 

The thermal conductivity of a unidirectionally reinfored composite 
in a direction parallel to the fibres X^ 1 is usually expressed in terms 
of the volume weighted conductivities of the constituents using an Ohm's 
law approach 

X 11 = X I* + (1-V_) X (8) 

c f f f m 

where X ^- 11 is the conductivity of the fibres 

X is the matrix conductivity 
m 

is the fibre volume fraction 

Applying the results for the three unidirectional composites should in 

principle permit the calculation of the thermal conductivity X ^ 11 of the 

fibres. However the significantly higher densities of the 1-D composites 

suggests that it would be inappropriate to use the measured matrix 

conductivity in equation 8 . Certainly the composite conductivity of the 

1-D fibre "F" composite (3.0V/ cm 1 K 1 at 300K) cannot be accounted for 

in terms of the measured matrix conductivity (0.85W cm 1 K 1 at 300K) 

and the reported conductivity of fibre "FP (0.6W cm 1 K 1 at 300K e 

The conductivity anisotropy ratios of the 1-D composites graphitized 

at 2750°C are approximately 11-12 at room temperature, decreasing to 

about 8 at 2300K. These are high values and irrespective of whether 

the bulk matrix properties are representative of the composite matrix, the 

implication must be that fibre anisotropy is also high. It has been 
( 22 ) 

reported that matrix conductivity is usally below both axial and 

transverse fibre conductivity. This is a surprising observation and are 

variance with our measured data. 

Various models have been put forward in an attempt to explain 

the transverse thermal conductivity of a unidirectional composite in 

terms of the thermal conductivities of the matrix and included fibres. 

(23) 

Some of these have been critically assessed by Pilling et al for 
fibre/epoxy composites. We adopt a modified Ohm's law approach, the 
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in-phase shear field analogy, using a model and analysis developed by 
(24) 

Springer and Tsai . This yields a generalised equation for the 1-L) 

X 

transverse thermal conductivity X 

J c 



where 0 = f(y), a function relating fibre width at any givoiy and X^ is 

the fibre transverse conductivity. The solution of Equation (9) is shown 

graphically in Figure 66, for a fibre of circular cross-section and square 

packing, as a function of fibre volume fraction (FVF) and X,. /X . Equation 

r m 

5 assumes good contact between fibre and matrix but Figures 34-37,43-57 have shown 
this assumption to be quite appropriate for the 1-D composite. 


The fibre volume fraction of the 1-D composites of fibre "A" is 
50%, that for the 1-D composite of fibre "F" is 52%. Using values for 
matrix conductivity calculated from diffusivity measurements and values 

determined for thermal conductivity in a direction perpendicular to fibre 

X 

orientation X we deduce, using data at 350K and 1000K. 


TABLE 3 


Temperature 

1 350K 

1 1000K 

spec type 

X 

X 

C A 

m 

X 

X f 

r /x 

m 

X 

X 

c /X 

m 

X 

X f 

m 

Fibre "F" no CVD 

0.33 

0.1 

0.38 

0.15 

Fibre "A" CVD 

0.12 

<0.01 

0.21 

0.01 

Fibre "A" no CVD 

0.23 

0.03 * 

0.26 

0.05 



There are indications that the ratio 'X is temperature dependent 
for all 1-D composites. 

If we now considered this data in more detail then three factors 
can affect the thermal conductivity of the 1-D composites. 


a) The density of matrix graphite in the 1-D composite is 
higher than that of the matrix graphite as supplied. 













b) The fibres within the composites will have been graphitised 
to a temperature of 2750°C. 

c) Bonding of matrix graphite to the fibre preform can confer 
some directionality on thermal properties for the composite. 

Matrix crystallites in the vicinity of non CVD'd type "F" fibres 

. (17) 

are oriented parallel to the fibre axis. 
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9. 3-D Composites 

With the present state of knowledge and lack of suitably relevant 
data the modelling of C.F.C.C. can only succeed if the problem of 
quantifying the influence of processing can be avoided. This is possible 
if the 3-D C.F.C.C. is considered to be a two phase material comprising 

a) Processed fibre yarns representing the exial reinforcement 

b) Isotropic matrix pockets filling the rest of the space in 
the unit cell geometry. 

9 ( i) C.F■C.C. 1 s A and C 

The method rests on the assumption that the yarns and matrix pockets 

are arranged either in series or in parallel and that the resulting thermal 

model is analogous to series and parallel connected electrical circuits. 

(26) 

Essentially, it extends the method used by Knappe and Martinez-Freire 
to three dimensions. 

Along the axis under consideration, the composite unit cell is sub¬ 
divided into four parallel conduction channels (see Figure 67;, jt is then 
required that: 

a) The temperature difference AT along the heat flow direction is constant. 

b) The total heat flow Q may be divided into four parts i.e. 

« = «L * «2 * 0, * \ (10) 

The thermal conductivity of each channel is calculated using an Ohm's 

law approach. Thus, for the Z axis channels indicated in Figure 67: 


giving 


2aX X, 
m lx 

(sX + gX ) 
Ax m 


where X^ x is the transverse conductivity of the X axis yarn. 
Similarly, 


2aX X 

Axly _ 

(sX + gX ) 

* ” JLy 


( 13 ) 
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2aA A 
m lj_ 


3Z (sA 1 gi ) 

m sy 


Lastly, 

X 4Z = X //Z 


(I'D 


(15) 


where A,, is the parallel conductivity of the Z axis yarn. 

The conductivities of the individual channels are then area-weighted 
and added to obtained the total conductivity of the unit cell. For the Z 
axis, the following result is obtained: 


X lZ^ 2b ~ h ^ r + X 2Z hr + X 3 Z^ 2 c- r)h + A 4 Z (2c-r)(2b-h) 

4bc 


(16) 


And for the X axis: 


X lX rS + X 2X^ 2c ~ r ^ S + X 3X^ 2c_r ^ g + X 4X rg 
4ac 


(17) 


2bA A 

where A^ = -. 


2hX XZj.Y 


[hA m t ( 2 b-h)^] 


2X 


[ hA lz + (2b-h)A xy ] 


2bA A 
Z m 


3X 


[hA + (2b-h)A ] 

j. z m 


* X 4X = X //X 


In the real composite, these results must represent only the upper 
bound of channel and unit cell conductivities. Account must also be taken 
of the effect of structural defects i.e. split yarns, broken interfaces and 
porosity, though the latter may already be adequately represented within 
the 1-D composite and bulk matrix data. 

For split yarns and broken interfaces perpendicular to heat flow, the 
lower bound to individual channel conductivity will be zero and from this the 
unit cell lower bound for the particular defect type can be found. A value 
between the bounds can then be determined from the frequency With which the 
defect appears in the composite structure. 

(ii) C.F.C.C. 'B' 
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and perpendicular to the weave of a 2-D CRCC, using a method due to 
Bruggeman. However, this approach uses the separately measured properties 
of the individual fibres and matrix. Our investigation shows that these 
properties are considerably affected by the CVD process and fibre bundle 
densification. 

Analyses of 2-D woven composites have agreed the weave structure has 

little influence on the conductivity perpendicular to the reinforcement 

plane'' ’ . Although the non orthogonal FWPF construction increases 

• (23 ) 

conduction path length, this is balanced by increased conductivity 

Figure 45 has shown that the X-Y weave is not significantly distorted 
by the Z axis yarns. It seems appropriate, therefore, to retain the 
orthogonal unit cell analysis for the Z axis conductivity of CFCC ’B'. 

The case of the parallel or in-plane (X axis) conductivity is more 

complex. The weave also results in increased in-plane conductivity paths, 

(23) ... 

but Kessler has estimated that the effect is negligible for the 

particular weave used in CFCC 'B * (eight-harness satin weave). Figures 45 

and 46 have shown the absence of clearly defined matrix cross-over pockets 

in the X-Y weave, the yarns tending to fill the whole volume available. 

The reduced incidence of the relatively high porosity matrix pockets 

may explain the higher density of CFCC 'B'. Given that the yarn axial 

conductivity is considerably greater than that of the bulk matrix, the 

orthogonal analysis must represent a lower bound to the X axis conductivity 

of this composite. * 

A simple upper bound can be determined from the assumption that the 

heat flow paths within the X-Y weave are confined to the X axis yarns. The 

low conductivity transverse paths through the Y axis yarns are then effectively 

'short circuited' . Using the orthogonal analysis, this may be implemented by 

expanding the X axis yarn into the volume normally occupied by the Y yarn- 

matrix conduction channel. Although this implies a reduction in yarn FVF, 

because of processing complications it is not clear that there must be a 

concomitant reduction in yarn conductivity such as would be implied by 







f 


Equation *1 • Equation 17 then becomes: 

A = ^ 2x (2c-r) t; + A 3X (2c-r) C + 2A^r 


4ac 










33 


10. Analysis of Data 

10.1 1-D C.F.C.C. Composites 

If we follow the approach suggested in 8.1 to deduce thermal 
conductivity values for the matrix graphite within the 1-D composite 
then a matrix density of 'v 1.9 gin cm 3 would be appropriate for the 1-D 
fibre "F" composite and 1.67 gm cm 3 for the fibre "A" composites. This 
would suggest porosity/tortuosity factors of 4.25 and 6.5 respectively 
and imply weighting factors of 2.1 and 1.4 for values of matrix 
conductivity to be used in equations 8 and 9. Applying this to equation 
8 yields room temperatures values of A^. 11 - 4.1 Wcm 3 K 1 for fibre "F" 
and A^ 11 =3.0W cm *K 1 for fibre "A" after graphitisation at 2750°C. 

It is not realistic to apply this approach to the CVD’d composite of fibre 
" A" since we have no value for matrix graphite graphitised to only 2300°C. 

Considering the perpendicular direction however and using these 
estimated values for matrix graphite conductivity of 1.8 Wcm' 1 K _1 for 

the fibre "F" composite and 1.2 Wcm” 1 K _1 for the fibre "A" composite 

J. 

yields values of A c /A m of 0.15 and 0.167 for the two composites respectively. 

The limiting value of A. /A < 0.01 would in turn suggest values for A 

r m m 

= 1.35 Wcm 1 K 1 in the fibre "F" composite and A m = 1.0 Wcm' 1 K _1 in 
the Fibre "A” composite. 

It must firstly be asked if these values are reasonable. We have 
concentrated on data at low temperatures because the only data available 
for comparison is low temperature dat’a (90-270 K) obtained for carbon 
fibre/epoxy resin composites * * . It is reasonable to suppose 

that the orientation of the "a" axis crystallites of a carbon fibre 

will be directly related to its thermal history. The data of Volga and 

(24) ( 28 ) . 
Pilling et al clearly shows this to be the case. Data by Volga 

show a dramatic increase in thermal conductivity of carbon fibres 

with graphitisation temperature. The fibre thermal conductivity A f n 

at 270 K, increases from 1.6 Wcm 1 K 1 after graphitisation for 1 hour 




at 2600°C to 3.0 Wcm 1 K 1 after 1 hour at 2800°C. Hence our derived 
values for are very much in agreement with these observations. 

Likewise Pilling et al have determined the transverse thermal 
conductivity of UTS carbon fibres manufactured from polyacrilonitrite 
precursor similar to fibre "A" to be 0.0G Wcrn 1 K 1 at 270K. This is 
again in accord with our conclusion, that the transverse conductivity 
of both fibre "A" and fibre "F" is at least one, and more nearly two 
orders of magnitude lower than the parallel conductivity. 

The one remaining uncertainty is why the ratio X^A m , using 

assumed matrix conductivity values, is less than the limiting value of 

0.2 predicted by the model for a F.V.F. circa 50%. Hitherto we have 

concentrated on using bulk properties and neglected the influence of 

(15) 

processing. Stover et al have shown that matrix crystallites in 
the vicinity of non C.V.D.'d fibre "F" composites are parallel oriented 
(P.O.G.) with respect to the fibre axis. Our electron microscopy 
observations neither prove nor disprove this observation. However if 
true this will result in an increase in effective fibre diameter and 
a concomitant reduction in the amount of randomly oriented matrix graphite 
in the 1-D composite which would adequately explain the low X c A ratio. 

10.2 3-D C.F.C.C. Composites 

A necessary prerequisite is to establish the compositional similarity 
between the 1-D composites and the 3-D yarns. In the case of the 1-D 
fibre a composite and composite "C" this is patently inappropriate since 

k 

the fibres in composite "C" are clearly of a different type. Nontheless 
since we have no accurate data for this composite we will proceed, for 
reasons that will be made clear later on the assumption that we can apply fibre 
"A''1-data to composite "C" 

For the general case the yarn fibre volume fraction is calculated 
from the cross sectional area, the number of filaments/yarn and the affective 
filament area. This yields a F.V.F. of 0.60 - 0.64 ^^which is higher than 
the F.V.F. of either 1-D composite. However the difference in processing will 
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result in a significant difference in properties thereby complicating any 

attempt to model the results. Since no experimental information is 

available to quantify the differences between CVD versus non CVD on fibre 

"F" composites, any attempt to determine the conductivity of a CVD'd yarn 

from the 1-D composite data can have doubtful merit only. However, since 

there are no 1-D composites with fibre "F", CVD and standard process 

( 6 ) 

available there is no alternative.Similar arguments do not now apply to 

composite "C". 

For fibre "F" 3-D composites it has been shown that 

a) Sheath zones primarily of P.O.G. can form around the fibres. 

b) Where filaments are closer together than about lg diameters the 
space may be entirely filled with C.V.D. graphite. 

c) For distances greater than this up to 7 fibre diameters the space 
is filled with transversely oriented pitch. 

d) At distances greater than approximately 7 diameters the orientation 
becomes random. 

It is further noted that where sheath zones predominate the matrix tends 
to align the crystallographic planes normal to the filaments. This presents 
a somewhat idealised concept of the unit cell of a 3-D composite as shown 
in Figure 65. Broken yarns are likely to be matrix filled. 

Our observations on the 1-D fibre composites indicate that the non 
CVD'd composites of fibre "A" and "F" show fairly dense composites with 
little evidence of preferred orientation in the matrix graphite. However 
the electron micrographs of the CVD'd composite of fibre "A" shows the formation 
of a distinct sheath around individual fibres and evidence for the formation 
of transversely oriented graphite in the region between fibres. The density 
of fibres is too high to determine whether the orientation of the matrix • 
graphite would randomise at increased fibre separation. 

It therefore seems logical that the change of orientation from 
P.O.G. to T.O.G. should significantly increase the transverse matrix 
conductivity. Again our arguments demonstrate that, although the 








thermal conductivity of matrix graphite is higher than that of the fibres 
in the transverse direction the latter will predominate. Tor example 
figure 66 shows that a 100a increase in matrix conductivity will result in 
a 35a increase in transverse conductivity of a composite with a similar 
F.V.F. A further complication concerns cracks, the density of which is 
higher in the 3-D composites than in the 1-D composite. The majority of 
these run parallel to the fibre axis. They should have little effect 
on the transverse thermal conductivity which is dominated by the matrix 
but have a significant affect on the axial thermal conductivity. 

On the basis of the foregoing discussion, the 1-D composite 
conductivity data is modified in the following way for use in the 3-D 
composite model. 

a) Axial conductivity reduced by a factor of 0.66 

b) Transverse conductivity increased by a factor of 2.0. 

No additional correction has been made for the difference in FVF 

noted earlier. The bulk matrix is used unmodified as its porosity is of a 

^ 0 \ 

similar order to that of the cross- over pockets v 
d) Modelling Results 


The appropriate unit cell constants are given in Table 4. 


Composite 

Unit Cell dimensions mm 

Type 

a 

b 

c 

g 

h 

r . 

s 

"A" 

0.42 

0.38 

0.38 

0.42 

0.32 

0.32 

0.42 

"B" 

0.127 

0.625 

0.625 

0.127 

0.71 

0.71 

0.127 

"C" * 

0.38 

0.43 

0.43 

0.38 

0.34 

0.34 

0.38 


* Used as supplied by AFML although this is known to be incorrect. 

Table 4 

Unite Cell Constants 


The predicted upper and lower bounds of C.F.C.C. "A" X axis conductiv 
are shown in Tables 5 and 6 together with the experimental data as defined 







by the least squares function. The rms errors are 27% and 11% respectively. 
Clearly the lower bound, representing the case of one broken yarn/yarn 
interface per unit cell, provides the best fit with the experimental data. 

It was apparent from the microstructural investigation (see Figures 38 and 39) 
that the lower bound assumptions more closely reflected the real material. 

The predicted upper and lower bounds of the CFCC "A" Z axis is shown 
in Tables 7 and 8. The rms errors are 12% and 9% respectively. The results 
are very similar to those of the X axis with the greatest error, in the case 
of the lower bound, occurring at the extremes of the temperature range. 

The predicted conductivity of CFCC "B" X axis is shown in Tables 
9 and 10. As discussed earlier, the upper and lower bounds of this material 
are different from those of CFCC "A". The rms errors are 10% and 20% respect¬ 
ively. 

As expected the lower bound, which represents the orthogonal unit 
cell geometry, consistently underestimates the X axis conductivity. The upper 
bound, which though exaggerating the weave influence, is in much closer 
agreement with the experimental data. 

The calculated upper and lower bounds of CFCC "B" Z axis 
conductivity are shown in Tables 11 and 12. The rms errors are 24% and 
14% respectively. The lower bound chosen is again that of one broken yarn/ 
yarn interface per unit cell and this seems reasonable . in view of the evidence 
provided by the microstructural investigation (See figures 45 and 46). 

For C.F.C.C. "C" we have used in the model data for fibre "A" 

> 

composites although it is clear that the fibres in composite "C" care not 
fibre "A" but more probably fibre "F". In Tables 13 and 14 the calculated 
upper and lower bounds are shown for the X axis. The rms errors are 6.6% 
and 22.8% respectively. The calculated conductivities for the Z axis are. 
presented in Tables 15 and 16 respectively and show r.m.s errors of 17.5% 
and 23.8% respectively. 
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Not surprisingly the modelling results for composite "C" are 
poorer than those for composites "A" and "B". The lower bound which is 
arguably the most reliable, in view of the microstructural evidence, 
and certainly gives a better prediction for composite "A" and the Z 
axis of composite "B", consistently underpredicts the thermal conductivity. 
Clearly higher conductivities particularly parallel to fibre orientation 
would be appropriate. Comparison of data for non CVD'd 1-D fibre "F" comp¬ 
osites (figures 24-27) and fibre "A" 1-D composites (figures 28-31) show 
the former to be clearly higher by some 50% at 350K decreasing to 25% 
at 1500K. If as we believe composite "C" is also a fibre "F" composite 
data based on figures 24-27 would be more appropriate in the model and lead 
to better predictions. 
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Conclusions 


Modelling of C.F.C.C. thermal properties can only succeed if 
accurately quantifiable information is available on the properties of 
the individual phases. In this case it has not been possible to rigorously 
test the present CFCC model because of the unavailability of any suitably 
processed 1-D composites. However, in spite of its basic simplicity 
and the limitations of the 1-D composite data, the predictions of the 
CFCC model have shown a reasonable agreement with the experimental data. 

For composites manufactured from fibre "F" material, with the 
exception of the X axis of composite "B" which is a special case CFCC 
conductivities are reflected more closely by the lower bound solutions. 

It is believed that the microstructural investigation has shown that these 
in turn reflect most accurately the structural state of the real composites. 
In the case of composite "C" microstructural evidence comparing this with 
1-D composites from fibre "A" leads to the conclusion that fibres are not 
type "A" but more likely fibre "F". This is reflected in modelling values 
which consistently under predict and it is suggested that better values 
might be obtained by using fibre "F" 1-D composite data. 
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e * am c 

LOWER pound 


TTMP£Rft7 jpe; \K) 

rxe. 04TA 

Ctl. DtlA 

fc-ER&OP 

30o.p 

1.783 

1.0H2 

-21.8 

350.0 

1.332 

1.055 

-20.7 

400.0 

1.283 

1.029 

-19.8 

450.0 

1.237 

1.003 

-18.9 

5oo.o 

1.194 

.978 

-18.1 

550.0 

1.153 

.95? 

-17.4 

600.0 

1.115 

.927 

-16.8 

650.0 

1.079 

.903 

-16.3 

700.0 

1.045 

.679 

-15.9 

750.0 

1.014 

.856 

-15.6 

fiOO.O 

• 9B4 

.833 

-15.4 

850.0 

.957 

.811 

-15.2 

900.0 

.931 

.790 

-15.2 

950.0 

.907 

.769 

-15.2 

1000.0 

• RB4 

.749 

-15.3 

1050.0 

.863 

.730 

-15.5 

1100.0 

.844 

.711 

-15.7 

115".0 

.825 

.694 

-16.0 

l?on.o 

.809 

.677 

-16.3 

1250.0 

.793 

.661 

-16.7 

1300.0 

.779 

.645 

-17.) 

1350.0 

.764 

.631 

-17.5 

1400.0 

.752 

.617 

-17.9 

1450.0 

.740 

.604 

-18.4 

1500.0 

.729 

.592 

-18.p 

1550.0 

.719 

• 5B1 

-19.3 

1600.0 

.710 

.570 

-19.7 

1650.0 

.701 

* .560 

-20.J 

1700.0 

.693 

.550 

-20.5 

1750.0 

.665 

.54? 

-20.9 

1500.0 

.679 

.533 

-21.3 

1850.0 

.87] 

.526 

-21.7 

1900.0 

.664 

.518 

-22.0 

1950.0 

.659 

.512 

-22.3 

?ooo.o 

.653 

.505 

-22.6 

2050.0 

.647 

.499 

-22.9 

2)00.0 

.64? 

.493 

-23.2 

2150.0 

.637 

.4B7 

-23.4 

2?00.0 

.632 

.4 6? 

-23.7 

2250.0 

.627 

.476 

-24.0 

2300.0 

.6?? 

.■471 

-24.3 

2350.0 

.617 

.465 

-24.7 

2400.0 

.612 

.459 

-25.1 

2450.0 

.606 

.452 

-25.6 

2500.0 

.603 

.446 

-26.1 


PMS ERonR= )9.9 


Table ID 








TnrKrtAL CO-<l'»iCl JVJTy K>Dn_l*J'. (>«■ CFCC 
« z amb 

UPP£R BOUND 


7rMp£RATUPE no 
300.0 
350.0 
400.0 
450.0 
500.0 
550.0 
600.0 
650.0 
700.0 
750.0 
BOO.O 
850.0 
900.0 
950.0 
1000.0 
1050.0 
1100.0 
1150.0 
1200.0 
1250.0 
1300.0 
1350.0 
1400.0 
1450.0 
1500.71 
1550.7!-- 
1600.0 
1650.0 
1700.0 
1750.0 
1800.0 
1850.0 
1900.0 
1950.o 
2000.0 
2050.0 
2100.0 
2150.0 
2200.0 
2250.0 
2300.0 
2350.0 
2400.0 

2450.0 

2500.0 


PMS EPROR= 24.1 


FXP. data cal. data 


.641 

.904 

.633 

.885 

.625 

.866 

• 616 

.847 

• 607 

• 828 

• 597 

.808 

.587 

.789 

.576 

.770 

.565 

.751 

.555 

.732 

.544 

.714 

.533 

.696 

.522 

.678 

.512 

.661 

.502 

.645 

.492 

.629 

.483 

.613 

.474 

• 59B 

.465 

.584 

.457 

.571 

.449 

.558 

.442 

.546 

• 436 

.534 

.430 

.523 

.425 

.513 

.420 

.503 

.416 

.494 

.413 

.486 

.410 

.478 

.407 

.471 

.406 

.465 

.404 

.458 

.404 

.453 

.403 

.447 

.403 

.442 

.404 

.437 

.405 

.433 

.406 

.428 

.407 

.424 

.409 

.419 

.410 

.415 

.412 

.410 

.413 

.405 

.414 

.400 

.415 

.394 


*£PPOR 

41.0 

39.7 
38.5 

37.4 

36.4 

35.4 * 

34.5 

33.6 

32.8 
32.0 

31.3 

30.6 

29.8 

29.1 
' 28.5 

27.8 

27.1 

26.4 

25.6 

24.9 

24.1 

23.3 

22.5 

21.7 

20.8 - 

19.8 

18.9 
17.R 
16.8 
15.7 

14.6 

13.4 

12.2 

10.9 

9.6 
8.3 
6.9 
5.5 

4.1 

2.7 

1.2 
-.4 

-1.9 

-3.5 

-5.2 


Table 11 





JMLRH&'- CO^PuO f v 1 Ty mode t_ I u“t yf Cf CC « 
B Z A*l< 

10W£R POUND 


TEMPERATURE l<) 

rxo. Data 

CAL. DATA 

if PROP 

300.0 

.64) 

.720 

12.4 

350.0 

.633 

.703 

11.0 

400.0 

.625 

.686 

9.7 

450.0 

• 616 

.668 

6.5 

500.0 

.607 

• 651 

7.4 

550.0 

.597 - 

.635 

6.3 

600.0 

.587 

.618 

5.4 

650.P 

.576 

.602 

4.5 

700.0 

.565 

.586 

3.6 

750.0 

.555 

.570 

2.8 

SOO.P 

.544 

.555 

?.l 

650.0 

.533 

.540 

1.3 

900.0 

.522 

.526 

.6 

950.0 

.512 

.512 

-.0 

1000.0 

.502 

.499 

~ -.7 

1050.0 

.492 

.485 

-1.4 

1100 .P 

.483 

.473 

-2.0 

1150.0 

.474 

.461 

-2.7 

1200.0 

.465 

.449 

-3.3 

3250.0 

.457 

.439 

-4.0 

1300.0 

.449 

.428 

-4.7 

1350.0 

.442 

• 41 B 

-5.4 

1400.0 

-436 

.409 

-6.1 

1450.0 

.430 

.400 

-6.9 

1500.-0 

.425 

.392 

-7.7 - 

1550.0 

.420 

.384 

-P.5 

1600.0 

• 41b 

.377 

-9.4 

1650.0 

.413 

.370 

-10.3 

1700.0 

.410 

.364 

-11.2 

1750.0 

.407 

.358 

-12.2 

1 BOO.0 

• 40b 

.352 

-13.2 

1850.0 

.404 

.347 

-14.? 

1900.0 

.404 

.342 

-15.3 

1950.0 

.403 

.337 

-16.4 

2000.0 

.403 

.333 

-17.5 

2050.0 

.404 

.329 

-IB.6 

2 )00.0 

.405 

.325 

-19.7 

2150.0 

• 40b 

.321 

-20.9 

2200.0 

.407 

.3)7 

-22.0 

2250.0 

.409 

• 3) 4 

-23.2 

• 2300.0 

.410 

.310 

-24.4 

2350.0 

.412 

.306 

-25.6 

2400.0 

.413 

.302 

-26.B 

2450.0 

.414 

.298 

-2B.0 

2500.0 

.415 

.294 

-29.3 


RMS FRRORs 13.6 


Table 12 





THERMAL 

conductivity 

MODELING or crcc 

C 


C X AXIS 



UPPER 

BOUND 


TEMPERATURE <K> 

EXP. DATA 

CAL. DATA 

TERROR 

350.0 

.624 

.658 

5.5 

400.0 

.616 

.647 

5.0 

450.0 

.608 

.636 

4.6 

500.0 

.600 

.625 

4.1 

550.0 

.591 

.613 

3.6 

600.0 

.582 

.600 

3.1 

650.0 

.573 

.588 

2.6 

700.0 

.563 

.575 

2.1 

750.0 

.554 

.562 

1.6 

800.0 

.544 

.550 

1.0 

850.0 

.534 

.537 

.5 

900.0 

.525 

.524 

-.1 

950.0 

.515 

.512 

-.6 

1000.0 

.506 

.500 

-1.2 

1050.0 

.497 

.488 

-1.7 

1100.0 

.488 

.477 

-2.3 

1150.0 

.479 

.466 

-2.8 

1200.0 

.471 

.455 

-3.4 

1250.0 

.463 

.445 

-3.9 

1300.0 

.455 

.435 

-4.4 

1350.0 

.448 

.426 

-4.9 

1400.0 

.441 

.417 

-5.3 

1450.0 

.434 

.409 

-5.8 

1500.0 

.428 

.401 

—6.2 

1550.0 

.422 

.394 

-6.6 

1600.0 

.417 

.388 

-7.0 

1650.0 

.412 

.381 

-7.3 

1700.0 

.407 

.376 

—7 • 6 

1750.0 

.403 

.371 

-7.9 

1800.0 

.399 

.366 

-8.2 

1850.0 

.396 

.362 

-8.4 

1900.0 

.393 

.359 

-8.6 

1950.0 

.390 

.356 

-8.8 

2000.0 

.388 

.353 

-9.0 

2050.0 

.386 

.350 

-9.2 

2100.0 

.384 

.348 

-9.3 

2150.0 

.383 

.346 

-9.4 

2200.0 

.381 

.345 

-9.6 

2250.0 

.380 

.343 

-9.7 

2300.0 

.380 

.342 

-9.9 

2350.0 

.379 

.341 

-10.1 

2400.0 

.378 

.339 

-10.3 

2450.0 

.377 

.338 

-10.5 

2500.0 

.377 

.336 

-10.7 


RMS ERROR* 6.6 


TABLE 13 







thermal conductivity modeling or crcc c 

C X AXIS 
LOVER BOUND 


TEMPERATURE(K) 

EXP. DATA 

CAL. data 

ftERROR 

350.0 

• 624 

.539 

-13.7 

400.0 

.616 

.530 

-14.1 

450.0 

.608 

.520 

-14.5 

500.0 

.600 

.511 

-14.9 

550.0 

.591 

.501 

-15.3 

600.0 

.582 

.491 

-15.7 

650.0 

.573 

.480 

-16.1 

700.0 

.563 

.470 

-16.5 

750.0 

.554 

.460 

-17.0 

800.0 

.544 

.450 

-17.4 

850.0 

.534 

.439 

-17.8 

900.0 

• 52S 

.429 

-18.2 

950.0 

.515 

.419 

-18.6 

1000.0 

.506 

.410 

-19.1 

1050.0 

.497 

.400 

-19.5 

1100.0 

.486 

.391 

-19.9 

1150.0 

.479 

.382 

-20.3 

1200.0 

.471 

.373 

-20.7 

1250.0 

.463 

.365 

-21.1 

noo.o 

.455 

.357 

-21.6 

1350.0 

.448 

.349 

-22.0 

1400.0 

.441 

.342 

-22.4 

1450.0 

.434 

.335 

-22.8 

1500.0 

.428 

.329 

-23.1 

1550.0 

.422 

.323 

-23.5 

1600.0 

.417 

.317 

-23.9 

1650.0 

.412 

.312 

-24.3 

1700.0 

.407 

.307 

-24.6 

1750.0 

.403 

.302 

-25.0 

1800.0 

.399 

.298 

-25.3 

1850.0 

.396 

.294 

-25.7 

1900.0 

.393 

.291 

-26.0 

1950.0 

.390 

.287 

-26.3 

2000.0 

.388 

.284 

-26.7 

2050.0 

.386 

.282 

-27.0 

2100.0 

.384 

.279 

-27.3 

2150.0 

.383 

.277 

-27.7 

2200.0 

.381 

.275 


2250.0 

.380 

.273 

- i C . _• 

2300.0 

.380 

.271 

-28. ! 

2350.0 

.379 

.269 

-29.1 

2400.0 

.378 

.267 

-29.5 

2450,0 

.377 

.265 

-29.9 

2500.0 

.377 

.262 

-30.4 



RMS ERROR* 22.B 


TABLE 14 











thermal conductivity modeling or cfcc c 

C Z AXIS 
UPPEP BOUND 


TEMPERATURE<K> 

EXP. DATA 

CAL. DATA 

ftERROR 

350.0 

.894 

• 825 

-7.7 

400.0 

.886 

.811 

-8.5 

450.0 

.877 

.797 

-9.2 

500.0 

.868 

.782 

-9.9 

550.0 

.857 

.767 

-10.5 

600.0 

.846 

.752 

-11.2 

650.0 

.834 

.736 

-11.7 

700.0 

.822 

.721 

-12.3 

750.0 

.810 

.706 

-12.8 

800.0 

.797 

.690 

-13.3 

850.0 

.784 

.675 

-13.8 

900.0 

.770 

.660 

-14.3 

95P.0 

.757 

• 646 

-14.7 

1000.0 

.743 

.631 

-15.1 

1050.0 

.730 

.617 

-15.5 

1100.0 

.717 

.603 

-15.9 

1150.0 

.704 

.590 

-16.2 

1200.0 

.691 

.577 

-16.5 

1250.0 

.679 

.564 

-16.8 

1300.0 

.666 

.552 

-17.1 

1350.0 

.655 

.541 

-17.4 

1400.0 

.643 

.530 

-17.6 

1450.0 

.632 

.519 

-17.9 

1500.0 

.622 

.509 

-18.1 

1550.0 

.612 

• 500 

-18.3 

1600.0 

.603 

.491 

-18.5 

1650.0 

.594 

.483 

-18.7 

1700.0 

.586 

.475 

-18.9 

1750.0 

.578 

.468 

-19.1 

1800.0 

.571 

.461 

-19.3 

1850.0 

.565 

.455 

-19.5 

1900.0 

.559 

.449 

-19.7 

1950.0 

.554 

.443 

-19.9 

2000.0 

.549 

.438 

-20.1 

2050.0 

.545 

.434 

-20.4 

2100.0 

.541 

.429 

-20.7 

2150.0 

.538 

.425 

-21.0 

2200.0 

.536 

.421 

-21.4 

2250.0 

.534 

.418 

-21.7 

2300.0 

.532 

.414 

-22.2 

2350.0 

.531 

.411 

-22.7 

2400.0 

.530 

.407 

-23.2 

2450.0 

.530 

.404 

-23.8 

2500.0 

.529 

• 400 

-24.5 


RMS ERROR* 17.5 


TABLE 15 








THERMAL 

conductivity modeling or crcc 

C l AXIS 

LOVER BOUND 

C 

TEMPERATURE IK) 

EXP. DATA 

CAL. DATA 

%ERROR 

350.0 

.894 

.763 

-14.6 

400.0 

.886 

.750 

-15.3 

450.0 

.877 

.737 

-16.0 

500.0 

.868 

.723 

-16.7 

550.0 

.857 

.709 

-17.3 

600.0 - 

.846 

.695 

-17.9 

650.0 

.834 

.681 

-18.4 

700.0 

.822 

.667 

-18.9 

750.0 

.810 

.653 

-19.4 

800.0 

.797 

.639 

-19.8 

850.0 

.784 

.625 

-20.3 

900.0 

.770 

.611 

-20.7 

950.0 

.757 

.598 

-21.0 

1000.0 

.743 

.584 

-21.4 

1050.0 

.730 

.571 

-21.7 

1100.0 

.717 

.559 

-22.1 

1150.0 

.704 

• 546 

-22.4 

1200.0 

.691 

.535 

-22.7 

1250.0 

.679 

.523 

-22.9 

1300.0 

.666 

.512 

-23.2 

1350.0 

.655 

.501 

-23.4 

1400.0 

.643 

.491 

-23.7 

1450.0 

.632 

.481 

-23.9 

1500.0 

.622 

.472 

-24.1 

1550.0 

.612 

.463 

-24.4 

1600.0 

.603 

.455 

-24.6 

1650.0 

.594 

.447 

-24.8 

1700.0 

.586 

.439 

-25.0 

1750.0 

578 

.432 

-25.2 

1800.0 

.571 

.426 

-25.5 

1850.0 

.565 

.419 

-25.7 

1900.0 

.559 

.414 

-26.0 

1950.0 

.554 

.408 

-26.3 

2000.0 

.549 

.403 

-26.6 

2050.0 

.545 

.398 

-26.9 

2100.0 

.541 

.394 

-27.3 

2150.0 

.538 

.389 

-27.7 

2200.0 

.536 

.385 

-28.1 

2250.0 

.534 

.381 

-28.6 

2300.0 

.532 

.377 

-29.1 

2350.0 

.531 

.373 

-29.7 

2400.0 

.530 

.369 

-30.3 

2450.0 

.530 

.366 

-31.0 

2500.0 

.529 

.362 

-31.7 


RMS ERROR* 23.8 


TABLE 16. 
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Figure 66 
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